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ABSTRACT 

We have supplemented existing spectra of Barnard's Loop with high accu- 
racy spectrophotometry of one new position. Cloudy photoionization models 
were calculated for a variety of ionization parameters and stellar temperatures 
and compared with the observations. After testing the procedure with recent 
observations of M43, we establish that Barnard's Loop is photoionized by four 
candidate ionizing stars, but agreement between the models and observations is 
only possible if Barnard's Loop is enhanced in heavy elements by about a factor 
of 1.4. 

Barnard's Loop is very similar in properties to the brightest components of 
the Orion- Eridanus Bubble and the Warm Ionized Medium (WIM). We are able 
to establish models that bound the range populated in low-ionization color-color 
diagrams (I([SII])/I(Ha) versus I([NII])/I(Ha)) using only a limited range of 
ionization parameters and stellar temperatures. 

Previously established variations in the relative abundance of heavy elements 
render uncertain the most common method of determining electron tempera- 
tures for components of the Orion-Eridanus Bubble and the WIM based on only 
the I([NII])/I(Ha) ratio, although we confirm that the lowest surface brightness 
components of the WIM are on average of higher electron temperature. 

The electron temperatures for a few high surface brightness WIM compo- 
nents determined by direct methods are comparable to those of classical bright 
H II regions. In contrast, the low surface brightness HII regions studied by the 
Wisconsin Ha Mapper are of lower temperatures than the classical bright HII 
regions. 

Subject headings: ISM:abundances-ISM:individual objects:Barnard's Loop- 
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Introduction 



The object commonly known as Barnard's Loop was di scovered photo graphically more 
than a century ago, originally by W. H. Pickering in 1889 ( jSheehanl 



E. E. 



19951 ) and then by 



Barnard! (118941 ) . Although a popular object for imaging by amateur astronomers, 



its low surface brightness has delayed the definitive observations necessary for developing 



a good model of its characteristics. A representative im age is shown in Figure 



f 

arc of about 14° facing east and at a distance of 440 pc (10' Dell fc Henney 



. It is an 



20081) is 110 pc 



in length. The northern boundary is at a Galactic Latitude of -13° and lines of galactic 
longitude are at a position angle of about 62°. These facts mean that its two sides (north 
and south) extend over a distance of about 60 pc and the conditions of the local ambient 
interstellar medium will have changed significantly. The north edge is about twice as bright 
as the south, indicating a higher density there. 

The structure of Barnard's Loop was the subject of a study of the brightness 
distribution of ultraviolet light measured i n a photograph obt ained during the Gemini 11 



manned space-flight. In this investigation 



O'Dell et al 



( 119671 ) derived the distribution of 



interstellar dust under the assumption that the ultraviolet continuum is due to scattered 
light originating in the OB stars in the Orion constellation Belt and Sword regions. It was 
found that the dust density increases approximately as the square of the dista nce from the 



center of Barnard's Loop. Under the assumption of a typical gas to dust ratio 



O'Dell et al 



( 119671 ) calculated that an ionization boundary should occur where one observes Barnard's 



loop to be brightest in Ha. It is expected that there woul d be neutral 



of this ionization boundary, which is consistent with older (IMenon 



lydrogen outside 



19581 ) and more recent 



1 Based in part on observations obtained at the Cerro Tololo Inter-American Observatory, 
which is operated by the Association of Universities for Research in Astronomy, Inc., under 
a Cooperative Agreement with the National Science Foundation. 
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( IHartmann fc Burton! 1 199 71 ) studies of HI 21-cm emissi on from this region. 



line splitting of abou t 24 km s 1 
the interpretation of |Q' Dell et al- 



i as bee n determined (IMadsen et al. 



Optical emission 



20061 ). consistent with 



( 119671 ) that Barnard's Loop is a result of the compression 
of ambient interstellar material pushed outward by radiation pressure force acting on the 
dust component, with the light originating from Sword and Belt region stars. 

Barnard's Loop is now understood to be only the bright eastern portion of a much 



l arger irregular she! 



([Reynolds fc Ogden 



of m aterial of 41° x 27° (380x220 pc) oriented east-west and thought 
19791 ) to be expanding at between 15 and 23 km s _1 . The entire object 



is called the Orion-Eridanus Bubble. The well delineated western part o 



Bubble shows the ion distribution of an ionization front (IMadsen et al. 



by neutral material (IHartmann fc Burton 



the Orion-Eridanus 



20061 ) surrounded 



19971 ) . The low velocity of expansion means that 



one would not expect mass- motion collisional excitation and ionization to be important, so 
that photoionization processes should dominate, but, the low velocities also make it d i fficult 



to discriminate between the radiation pressure driven model invoked by 



O'Dell et al 



(1967) 



and a mass-loaded shell o 



to explain the expansion (IMadsen et al. 



a supernova remna nt, the interpretation most commonly applied 



20061 ). 



The conditions of the gas i n the Orion-Eridanu s Bubble, as indicated by emission 



line ratios, appears to resemble (IMadsen et al 



20061 ) those of the very low-density, hot 



component of the interstellar medium called the WIM (Warm Ionized Medium) and contrast 
with those in the bright Galactic H II regions. Therefore, explanation of the conditions 
within the Orion-Eridanus Bubble may lead to understand ing the heating pro cesses within 



the WIM, which are presently not adequately understood (IHaffher et al. 



2009J). 



In this paper we will (§ 2) summarize the results from other studies of Barnard's Loop 
and present new spectrophotometric observations, then use these observations to test (§ 
3) models for the photoionization of the Orion-Eridanus Bubble, demonstrating that the 
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nebular lines, except for [0 III], can be explained by photoionization by the best four 
candidate ionizing stars only if Barnard's Loop is enhanced in heavy elements by about a 
factor of 1.4. In § 4 we argue that the low surface brightness ubiquitous component of the 
interstellar medium (ISM) commonly known as the WIM is subject to the same physical 
processes as the Barnard's Loop, that the electron temperatures in WIM components 
are similar to those of the classical bright H II regions, and that the common method of 
determining exact electron temperatures (T e ) of WIM components is rendered uncertain 
because of variations in heavy element abundances comparable to those seen in well studied 
H II regions. 

2. Observations 

Due to its low surface brightness Barnard's Loop has not been extensively studied with 
adequate spectral resolution. It is a popular object for imaging with wide-field cameras but 
usually with a filter that passes both the Ha emission line at 6563 A and the adjacent [N II] 
doublet at 6548 A and 6583 A, making it impossible to develop a quantitative analysis. The 
fainter portions of the Orion-Eridanus Bubble have been even more difficult to study. In 
this section we will summarize previous observations of the region and present new material 
on Barnard's Loop. 

2.1. Spectrophotometry of Barnard's Loop 

2.1.1. Previous Observations 

There have been two earlier studies that are particularly useful. The older study was 
at low spectral and high spatial resolution and the more recent at quite high spectral 
resolution but low spatial resolution. 
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A photoelectric scanner study of several faint regions of the sky (jPeimbert et al.lll975l ) 
included useful results for one position in Barnard's Loop (shown in Figure 1) with two 
5'.'2 x 77? 6 slits separated by 168". The spectral resolution is not stated but was adequate 
for separating the Ha and 6583 A lines. Although of apparently low signal to noise ratio 
(only the brightest lines were detected), this study extended to the bright [O II] doublet 
at 3727 A. The average surface brightness of these two samples was 1.1 x 10 7 photons 
cm -2 s _1 steradian -1 in Ha or 140 Rayleighs (the latter units are often employed in studies 
of extended emission and are 47r/10 6 times the surface brightness expressed in photon rate). 
The line flux ratios are presented in Table 1. 

A scanning Fabry-Perot spect rophotometer ( the W isconsin Ha Mapper, known as 



WHAM) study was conducted by iMadsen et al.l (120061 ) with a 1° diameter field and a 



spectral resolution of 12 km s , measuring the [O III] 5007 A, Hel 5876 A, Ha 6563 A, 
[N II] 6583 A, and [S II] 6716 A lines. The Ha surface brightness of position 1 was 86.3±0.1 
Rayleighs and for position 2 it was 228.1±0.2 Rayleighs. The relative fluxes are given in 
Table 1. 



2.1.2. New Spectrophotometry 

New spectroscopic observations were made with the Cerro Tololo Interamerican 
Observatory 1.5 m-telescope operated in cooperation with the SMARTS consortium. The 
instrument was the Boiler and Chivens spectrograph. Observations were made on 2008 
November 24 with Grating G58 using the Loral IK CCD detector. A GG395 glass filter 
was used to prevent second order flux from contaminating the first order flux that was 
targeted. One pixel projected 1'.'3 along the 429" long slit, while the slit width was 2'.'15. 
The measured full width at half maximum intensity of the nebular lines were 6 A. The 
CTIO spectrophotometric standard star Feige 15 was observed nine times and the results 
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were used to calibrate the nebular observations into energy units. The positions of the slit 
setting is shown in Figure 1 in addition to the position used for determining the background 
sky brightness. The sample was 6?07 from the Trapezium stars that contain 8 1 Ori C, the 
earliest spectral type star in the Orion Belt and Sword region. 

In order to facilitate sky subtraction and cosmic ray cleaning double exposures were 
made in the pattern of two 1800 s Barnard's Loop exposures followed by two 900 s sky 
exposures, a pattern repeated three times. Since it was considered most important to 
produce good sky subtractions near the Ha line, we scaled the sky observations to null out 
the OH band sky signal close to Ha prior to making the sky subtraction. This meant that 
other strong night sky lines that varied in brightness relative to the night sky OH emission 
were inadequately subtracted from the nebular spectra. 

Data reduction was done using standard IRAF procedures^. The results of these steps 
were calibrated spectra expressed in ergs cm -2 s _1 pixel -1 that were then converted to 
surface brightness units and averaged over the entire length of the entrance slit. Only the 
H 7 4340 A, H/3 4861 A, Ha 6563 A, [N II] 6583 A, [N II] 6548 A, [S II] 6716 A, and [S II] 
6731 A lines were bright enough to be measured using the task "splot", which required 
de-blending the lines near Ha and the [S II] doublet using task "deblend". The results 
from two spectra of the central Orion Nebula obtained during the same observing run were 
used to compare the derived surface brightness in the H/3 line with those obtained from the 



spectrophotometric study of 



Baldwin et al. 



telescope WFPC2 emission line filters 



19911) and the calibration of th e Hubble Space 



O'Dell fc Doil (119991 ). which use the 



Baldwin et al. 



(Il99ll ) results as a standard. There was good agreement. The surface brightness in H/3 was 



2 IRAF is distributed by the National Optical Astronomy Observatories, which is oper- 
ated by the Association of Universities for Research in Astronomy, Inc. under cooperative 
agreement with the National Science foundation. 
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5.3 x 10 6 photons cm -2 s _1 steradian -1 or 66.6 Rayleighs. The final averaged spectrum 
is shown in Figure 2. The negative values at 5577 A, 6300 A, and 6363 A reflect the 
uncertainties in subtracting the sky background as these [O I] sky lines varied considerably 
in surface brightness and the sky signal at 5577 A was comparable to the nebula's signal 
at Ha 6563 A. An upper limit for the nebula's [O III] 5007 A and He I 5876 A lines was 
determined by scaling the B./3 line and inserting it at the locations of those lines. The flux 
ratio of the 6716 A and 6731 A lines was within 1% of the theoretical low density limit of 
1.50. The relative line intensities, normalized to H/3 are given in Table 1. 



2.2. Earlier Observations of the Orion-Eridanus Bubble 

The 



Madsenetal 



(120061 ) study also included high signal to noise ratio fixed-pointing 
samples at four additional positions to the east of the center of OB stars in the Belt and 
Sword regions of the nebula, extending out to a distance of 24?0. They also obtained 
seven closely spaced samples crossing an ionization front at the southwest boundary of 
the Orion-Eridanus Bubble and lower quality data in regions further from the Galactic 
Plane and outside but near the Orion-Eridanus Bubble, thus sampling the more general 
WIM. When they compared their results with a similar set of measurements of large 
and low density H II regions they noted that the [S II] /Ha and [N II] /Ha flux ratios 
for the Orion-Eridanus Bubble varied together in the same way as in the WIM and they 
concluded that the variations in flux ratios were due to variations in T c , with the WIM 
being significantly warmer than the Orion-Eridanus Bubble material. They established 
that the pattern of the line ratios was quite different for the H II regions they had sampled 
and the Orion-Eridanus Bubble and WIM. The difference was in the sense that the [S II] to 
Ha ratio was much weaker in the H II regions at the same value of the [N II] to Ha ratios. 
Moreover, they noted that the [S II] 6716 A over [N II] 6583 A flux ratio was both larger 



than in their H II regions and increased only slightly with surface brightness in Ha, with an 
average value near 0.8. The value of this ratio was about 0.4 in their low surface brightness 
H II regions. 



3. Modeling 

It is important to determine the physical conditions in the Barnard's Loop and the 
Orion-Eridanus Bubble in order to understand both the nature of these objects and their 
origin. As we will see, there is a fundamental problem when trying to explain the line ratios 
(which reflect the conditions of ionization and excitation) with the mechanism of direct 
stellar photoionization by the hottest star in the region, but, we have been able to develop 
a model consistent with all the available observational material. 

There are two methods of testing the source of illumination of Barnard's Loop. The 
first test is to determine if the total stellar luminosity in the LyC that is derived from 
the surface brightness is consistent with the available flux from the enclosed stars and is 
considered in § 3.2. The second test is to determine if the observed emission-ine ratios can 
be explained by photoionization caused by the enclosed stars and is considered in § 3.4 for 
the well studied object M 43 as a test object and then in § 3.5 we apply this method to the 
Barnard's Loop. 



3.1. Basic Constraints Identified from the Observations 



The observations of the emission distribution, the surface brightness, and the known 
distance allow us to determine some basic properties of Barnard' s Loop. Fortunately 



the interstellar exti nction in the region is known to be negligible (IPeimbert et al 



1975 



Madsenetal 



20061 ). a conclusion consistent with the fact that our new Ha/H/3 flux ratio 
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of 2.90 is that predicted for an unreddened gas of about 9000 K. 



Baldwin et al. 



( 1199 ll ) demonstrated 



In their fundamental study of the Orion Nebula, 
that the surface brightness in the B./3 recombination line along a line of sight from the 
observer to the photoionizing star would be proportional to the flux of ionizing photons 
0(H) in units of photons cm -2 s _1 reaching a thin slab of gas. This approximation is 
probably valid for Barnard's Loop because the low inner-region densities would not be 
expected to absorb LyC photons. If one views the emitting layer face-on, then the relation 
would be </>(H)=47r(o!B/ttH^)S(H^) where the surface brightness S(H/3) is expressed in 
photons cm -2 s" 1 steradian -1 and a-Q and respectively the total recombination 

coefficient for hydrogen and the effective recombination coefficient for the H/5 line. The 
ratio of the recombination coefficients is not very sensitive to the assumed T e . 

In the case of Barnard's Loop we are observing a relatively thin shell seen nearly 
edge-on, so that the surface brightness will be enhanced by a factor of g above that for 
a slab viewed face on. We call the quantity the limb-brightening correction. In this thin 
shell model the surface brightness distribution would start at the maximum radius R max , 
rapidly rise to a peak value at angle corresponding to the inner radius R pea k, then slowly 
decrease to the value that would apply along the line of sight through the ionizing star. The 
enhancement of the surface brightness will be g=2( Rmax'Rpeak^VlR-max - Rpeak)- Using 
the dimensions for the line from 8 1 Ori C through our CTIO-2008 observations (R max =7?07 
and R pea k=6?01), the geometric correction factor is g=7.03, that is, the observed peak 
surface brightness will be enhanced by a factor of 7.03. 

One can derive the flux of ionizing photons (0(H)) from S(H/3) after the former has 
been corrected for limb brightening. Using our observed surface brightness (5.3 x 10 6 
photons cm -2 s -1 steradian -1 ), the geometric correction factor of 7.03 yields a flux of 
ionizing photons 0(H) =8.1 x 10 7 photons cm -2 s" 1 . It is hard to ascribe a probable error 
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to this number. The observational uncertainty is less than 10%. The uncertainties in the 
derived results that are due to applying a simple thin-shell model to Barnard's Loop, which 
has some internal structure, are probably significantly larger. As we will see in the next 
section, this flux and its wavelength distribution determine the expected photoionization 
structure and observed line ratios within Barnard's Loop. Having calculated the ionizing 
flux in a given sample, one can calculate the total luminosity in LyC photons (Q(H) with 
the units photons s _1 ) for the ionizing star(s) by the relation Q(H)=47r R 2 0(H). Adopting 
the average of the maximum and peak surface brightness radii of 1.56 x 10 20 cm yields 
Q(H)=2.5 x 10 49 photons s" 1 . 

This Q( H) value is larger than the value for the Orion Nebula of 7.8 x 10 48 photons 



s 1 found by 



Peimbert et al. 



(119751 ) from ear 



for the Orion Nebula of 1.1 x 10 49 found by 



y radio continuum o 



Dservations and the value 



van der Werf fc Gosd (119891 ) from their VLA 



study. The Orion Nebula Q(H) values are in approximate agreement with those expected 
from dominant photoionization of that object being by 9 l Ori C and # 2 Ori A, where the 
expected Q(H) for the s tars are 6 x 1 48 ph otons s _1 and 1.5 x 10 48 photon s" 1 respectively, 



using the calibration of 



Heap et al. 



( 120061 ). The fact that the predicted stellar value of 
Q(H) is smaller than that derived for the nebula probably indicates the uncertainty in the 
calibration or that the true spectral type is slightly earlier. 

One can also estimate the electron density from a knowledge of the surface brightness, 
the emissivity in the H/3 line and the geometry. For the constant density shell model the 
electron density (n e ) will be n 2 =47r S(H/3)/[a^ g(R m ax-R P eak)], assuming that all the free 
electrons arise from the photoionization of hydrogen, a very good approximation in this 
low-ionizationization region where helium is neutral (as indicated by the very weak or 
abse nt Hel 5876 A line). Us i ng th e previous values and adopting a^=3.63 x 10 -14 cm 3 s" 1 



from 



Osterbrock fc Ferlandl ( 120061 ) yields n e =3.2 cm 3 for the region of Barnard's Loop that 
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Madsen et al. 



we have observed and 0.7 cm 3 for the 
are similar to the density n e =2.0 cm -3 derived by 



(2006) western arc. These numbers 



Heilesetal. 



(120001 ) from combined optical 



and radio data over large samples in Barnard's Loop. 



3.2. Comparison of the Derived LyC Luminosity with that of the Enclosed 

Stars 



A natural first step in explaining the observed emission line properties is to investigate 
whether the line ratios can be expl ained by photoionizati on from a single dominant star, an 



approach adopted qualitatively by 



Peimbert et al. 



JlSZg). Wewi 



except that we shall use the sp ectral types summ arized in 



1 follow their discussion 



values from the recent study by 



Heap et al 



Goudisl Jl982j) and use the Q(H) 



(120061 ). There are many candidate ionizing stars 
within Barnard's Loop and its extension into the Orion-Eridanus Bubble, the brightest 
six in the LyC are the brightest Trapezium 06 V star 9 1 Ori C with an expected LyC 
luminosity of 6 x 10 48 photons s _1 , the cooler (spectral type 09.5 II) S Ori with an expected 
LyC luminosity of 5.6 x 10 48 photon s _1 , £ Ori (09.5 lb) at 5.6 x 10 48 photon s" 1 , and i 
Ori (spectral type 09 III) at 6.6 x 10 48 photon s _1 , plus # 2 Ori A and a Ori (both 09.5 V 



and 1.5 x 10 photon s ). It is unclear if e Ori belongs in this li st. With a spectra. 



Heap et al. 



( I2006h . 



Vaccaetal 



(119961 ) did 



type 



of BO la it is cooler than the stars studied by 
include one star (HD 37128) of this spectral type in their study and found it to be be 
6000 K cooler than an 09.5 la spectral type star (HD 30614) and would therefore have a 
much lower LyC luminosity. Since £ Ori has a spectral type of 09.5 lb, we conclude that e 
Ori's LyC luminosity is much lower than 5.6 x 10 48 photon s _1 and we will not include it 
in our tally. The location of these stars are shown in Figure 1. Although there are many 
additional stars of later spectral type and less luminous stars are distributed throughout the 
inner Orion constellation, their contribution to photoionization of Barnard's Loop must be 



-14- 



minimal. For a summary of the properties of these stars see Table 1.1. IV of iGoudisi (Il982l ). 



The H/3 surface brightness yields an ionizing flux of 0(H)=8.1 x 10 7 photons cm -2 
s _1 that then leads to a derived ionizing luminosity of Q(H)=2.5 x 10 49 photons s -1 . 
Photoionization of the large scale Barnard's Loop may arise from multiple stars. It is 
unlikely that two candidate stars that are located within the Orion Nebula (8 1 Ori C 
and 9 2 0i\ A) are im portant contributors to ionizing Barnard's Loop because it is well 



known (jO'Dell 



20011 ) that the nebula is optically thick to LyC radiation in all directions 



except possibly to the south west, where the fo reground Veil thins and X-ray emission from 



million degree gas is found ( jGiidel et al. 



20081 ) .The four remaining candidates for causing 



photoionization would have a total luminosity of Q(H)=1.9 x 10 49 photons s _1 , in reasonable 
agreement with the LyC luminosity derived from the FL3 surface brightness. The effective 
temperatures of these stars range from 30000 K to 32000 K. 



3.3. Method of Calculating Photoionization Models 



We have determined the properties of the radiation field illuminating Barnard's Loop 
by comparing the predictions for emission line ratios from stellar atmosphere models of 
various effective temperatures (T star ). W e use the deve l opme nt version of the spectral 



simulation code Cloudy, last reviewed by 



Ferland et al. 



(119981 ) and noted in appendix A. 



We present here several large grids of model calculations that, in some cases, varied both 
T s tar and the flux of hydrogen-ionizing photons, and in others, only T sta r- Later we also 
calculate and present models including differences in the abundance ratio ( Z/H) of heavy 



eleme nts to hydrogen. In many ways our calculations are similar to those of 



Sembach et al 



( 120001 ) but we use a version of Cloudy with up-to-date atomic coefficients, more recent 



stellar atmosphere models, and more recent calibrations of the spectral type-T sta r relation. 
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The other key parameters in the calculations are the abundances of the heavy elements 
that dominate the cooling of the ionized gas. We adopted N/H=6.5 x 10~ 5 , 0/H=4.3 x 10~ 4 , 
and S/H=1.4 x 10 -5 for these as a point of reference, using these for the M 43 calculations, 
where these abundances are known to apply, and scaled the abundances for the Barnard's 
Loop and WIM calculations, where there are no direct determinations of t he abundance s. 



The O/H and N/H ratios are a verages of the similar interstellar medium (jjenkinsl 120091 ) 



study of 



and M 4 2 (ISimpson et al 



Daflon et al 



2004 ) results and the S/H ratio of 1.4 x 10 5 was taken from the 



(120091 ). who measured two stages of ionized sulphur in ten B stars in 
the Orion Association. We also adopted a ratio of extinction to reddening of R=3.1, the 
standard dust to gas r atio, and included PAH's, and the TLUSTY stellar atmospheres of 



Lanz fc Hubenyl (120031 ). Density is a less important factor, but we have used the appropriate 



values for the two specific objects being modeled. 

We must note that the S/H ratio derived f rom H II regions can be significantly s malle r 



than the solar abundance of S/H =1.4 x 10 5 (lAsplund et al. 



20061 ). 



Simpson et al 



( 120041 ) 



give S/H=7.0 x 10 6 fro m IR lines in t he Or ion Nebula, a value consistent with previous 



optical studies. However 
the solar value, and the 



Daflon et al. 



Esteban et al. 



(2009) find that Orion B stars have S/H close to 
(120041 ) study of optical emission lines find relatively 
high values for the Orion Nebula ranging from 1.1 — 1.7 x 10~ 5 . Jenkins (2009) finds S/H 
depleted below the solar value in the ISM. These S/H values range over a factor of two, and 
Jenkin's (2009) discovery of a depletion pattern in the ISM means that S/H may depend 
on the location in the Galaxy. This affects Figure 5. The regions where [S II] lines form 
are mainly cooled by lines of [S II], [N II], and [O II]. We have adopted S/H = 1.4 x 10~ 5 . 
Had we used a lower value the [S II] lines would be weaker, but by l ess than the change in 



the S /H ratio because of the thermostat effect in a photoionized gas (jOsterbrock fc Ferland 



20061 ). If the [S II] cooling is diminished by lowering S/H, the gas will grow hotter, making 



lines of [O II] and [N II] stronger. The effect would be to shift the curves in Fig 5 down 
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and to the right. A factor of two decrease in S/H, which would be consistent with the older 
emission line studies , would shift the log 0(H) = 6.7 curve to roughly the position of the 



(IMadsen et al. 



by 



Simpson et al 



2006 1) dat a on this plot. A similar large range in N/H abundances was found 



19951 ) in their study of galactic H II regions, something that becomes 
important in our discussion of the indirect method of determining T e in WIM components. 
These known abundance variations play an important role in our interpretation of line 
ratios of both the Barnard's Loop and the WIM, as discussed in § 3.5 an d § 3.6. Our 



abund ance ratios most closely resemble the B star abundances adopted by 



Sembach et al 



(120001 ) except that we assume a higher overall abundance. 



The interstellar medium on the galaxy-wide scale is likely to be inhomogeneous and 



porous that is, certain 
by nearby dense clouds 



regions 



may 



De relatively clear of material, and others affected 



Spitzerl (119781 ). Quite sophisti cated mode 



Wood et al. 



s have been developed 



fl2005h : 



Haffner et al- 



to tak e this geometry into account (see, for example 

( 2009L . Although Cloudy is capab le of simulating complex geometries (the code Cloudy 



3D, see iMorisset fc Stasinskal (120081 ) and http://sites.google.com/site/cloudy3d/), here we 
use it to simulate a shell of gas that is symmetric about the central star or star cluster. 
The temperature, ionization, and emission in many thousands of lines are determined as 
a function of distance away from the central ionizing stars. Line and continuum optical 
de pths are taken into account so that changes in the ionizing radiation field (Chapter 



3, (jOsterbrock fc Ferland 



20061 )) are reproduced and the geometry is an inhomogeneous 



mix of temperature, ionization, and emission. In this mode a shell of gas, which appears 
to describe both the Barnard Loop and M 43, is closely mimicked. The quantities we 
report correspond to the observed spectrum along a tangent occurring at various impact 
parameters through the shell away from the star cluster. In this way we closely mimic the 
observations, which place a spectrometer slit at various positions. 



-17- 



There are at present no reliable dielectronic recombination rate coefficients for 
recombination from S +2 to S +1 . We reply on empirical estimates of the rate co efficient 
based on photoionization modeling of astronomical observations, as was done by 



Metal 



(119911 ) . Based on such models we would judge our current estimates to be uncertain by 



about 30%. This introduces roughly a 20% uncertainty in the balance between [S II] and 
[S III]. This produces an uncertainty that enters as an unknown scale factor that affects all 
models. This represents a basic uncertainty that affects all observations by this scale factor. 
It may shift the model predictions by an unknown systematic amount, but will not produce 
object to object fluctuations in the spectrum. 



3.4. Derivation of the Ionizing Star Properties from Emission Line Ratios in 

M 43 

To demonstrate the validity of the method employed in this paper to Barnard's 
Loop, we first made calculations that simulate the H II region M 43 (NGC 1982), which 
is dominated by a single cool star. We then used the results of that study to guide the 
calculations of the conditions for Barnard's Loop. 

M 43 lies to the immediate northeast of the much brighter Orion Nebula (M 42, 



O'Dell & HarriJ ( 


2010) 


(Schild & Chaffee 


1971 



Penston et al. 



1975 



Simon-Diaz et al. 



20081 ). The total Lyman 
continuum luminosity of the nebula is about 4.7 ± 2.0 x 10 47 photon s s -1 as determined 



from multiple radio wavelength studies summarized by iGoudisi (1l982l ). The apparent radius 
is 128" and the distance is 440 pc, giving a true radius is 0.27 pc. If the radius is adopted as 
the distance between the central star and a blister-type ionization front as applies to M 42, 
then the corresponding 0(H) value is 5.3 x 10 10 photons cm -2 s -1 . We have adopted an 
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electron density in M 43 of 520 cm" 3 flO'Dell fc Harris! l2010h . We do not use the ionizing 
luminosi ty expected from the spectral type of NU Or i because it is at least a triple star 



system (IMorrell &: Levato 



1991 



constant strong X-ray flux ( 



Preibisch et al 



Stelzer et al. 



19991 ) and is a known windy star with a 



20051 ). Therefore, we take advantage of knowing 



the total luminosity of the nebula and leave the temperature of the ionizing star as a 
variable. 

In this case we assume that the abundance of elements is that of Section 3.3 but 
that we don't know the exact geometry (the apparent radius need not be the same as 
the separation between the central star and the ionization front), so that we have varied 
the adopted $ values for a series of models for stars of varying T star (27542 K through 
40040 K in steps of 0.0325 dex), all with stellar gravities corresponding to main sequence 
stars. The results are shown in Figure 3. To these color- color diagrams we have added the 



observed line ratios for the six samples in M 43 from the 10' Dell fc Harris! (120101 ) study that 
gave I([N II] 6583 A)/I(Ha)=0.41±0.02, I([S II] 6716 A+6731 A)/I(H«)=0.18±0.03, and 
I([0 III] 5007 A)/I(H/3)=0.26±0.11. 

Examination of the low-ionization color-color diagram (Figure 3-left) shows that 
the best agreement of the emission line ratios with the predictions of the model stellar 
atmospheres is for a star with T star between 29700 K and 34500 K and log $~11.0. These 
temperatures bra cket the temperature of ~30000 K expected for a spectral type B0.5V 



( IHeap et al 



2006! ) and the best log $ value indicates that the distance between the ionizing 
star and the ionization front is 0.7 times the app arent radius of the nebula. A B0. 5 V 

20061 ) of 8 x 10 47 



star would be expected to have a LyC luminosity ( lOsterbrock fc Ferland 
photons s _1 , within a factor of two of the luminosity derived from the radio continuum. 

Examination of the high- ionization color-color diagram (Figure 3-right) shows that the 
I([0 III] 5007 A)/I(H/3) observations are best fit by a much higher T star . The cause of this 
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is probably not in the models, rather that the line ratio is co ntaminated by light or iginating 



in the nearby Huygens Region, where this ratio is about 3 (jO'Dell fc Harris 



20101 ). and it 



would only take about 20 % contribution to I([0 III] 5007 A) to account for the discrepancy 



O'Dell fc Harris! ( 120101 ) argue that essentially all of the ionization of M 43 arises from NU 
Ori, but their analysis does not preclude this level of contamination by scattered light. At 
similar distances within M 42 the scattered light from the Huygens Region is about this 
level. 



va 
in 



We show in Figure 4 the results of calculations for predictions of T e for the most likely 
ue of log <fr determined from Figure 3-left. We have superimposed the average T e derived 



O'Dell k Harrisl (120 10f ) from the nebular to auroral transitions of [N II] (7950±720 K) 



and [S II] (7260±720 K). The uncertainties in the temperatures (because of the weakness 
of the auroral transitions) is too large to narrowly define the stellar models, but we see that 
again the T e taken together agree with T star derived from the spectral type. 



3.5. Derivation of a Model for Barnard's Loop from the Ratios of the 

Strongest Emission Lines 

We then performed a similar set of calculations for Barnard's Loop, using the 
parameters relevant for that object. Fortunately, the constraining factors are better known 
than for M 43. In this case we know that all of the candidate ionizing stars (5 Ori-09.5 II), 
C Ori-09.5 lb, i Ori-09 III, a Ori 09.5 V) fall in a narrow range of T star near 31000 K. The 
side-on view of Barnard's Loop means that we know the separation of the ionizing stars 
and the ionization front reasonably well and in § 3.1 derive an ionizing flux of $=8.1 x 10 7 
photons s _1 . We also know that the density is about 3.2 cm -3 for our sample. Our 
calculations were made using T star =31000 K and a stellar gravity appropriate for a giant 
star. We adopted the previously assumed properties of interstellar grains and PAH and 
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initially used the same abundaces as in our calculations of M 43 models. 

In these new calculations there was a significant disagreement between the predicted 
model and the observations in both the low-ionization color-color (Figure 5) and the 
high-ionization color-color (Figure 6) diagrams. This disagreement was present whether 
one uses all of the available observations or only the ratios from this study. We tried to 
fit models with different values of $ and T star as low as 27500 K, but none of these agreed 
with the observations. The locus of points for various values of log $ for the closest low 
T sta r=27500 K models is also shown in Figure 5 . The small discrepancy at 27500 K 
becomes larger at higher values of T star since plots for T star progress to the right in the 
diagram, with the T star matching the likely ionizing stars shown as an open square around 
the a filled circle (log $=7.9). 

We then explored the results of varying the only parameter in our set of assumptions 
that is not directly determined for the Barnard's Loop, that is, the relative abundances of 
the heavy elements (Z/H). We scaled all of the heavy elements together, thereby ignoring 
subtleties in things like the N/O ratio changes that can be expected from stellar evolution 
models. In this case we see that one obtains excellent agreement between the observations 
and the models with a log (Z/H) abundance enhancement of between 0.1 and 0.2, i.e. a 
heavy element abundance enhancement of about a factor of 1.4. 

The more diagnostically useful figure is the low-ionization color-color ratio because of 
its insensitivity to contamination by scattered light and Figure 5 will be used for discussion. 
It is important to understand that the locus of points tracked by the models of various 
abundances is a result of both the number of elements of that species available for emitting 
the emission line, but also T c (the emissivity of the collisionally excited [S II] and [N II] 
lines increasing with T e and the emissivity of the recombination Ha decreasing with T c ). A 
track from lowest Z/H to highest Z/H models is a monotonic progression from high to low 
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Comparison of our models with varying abundan c es with the three set s of observations 



of Barnard's Loop (this study, (IPeimbert et al.lll975l ). (IMadsen et al.ll2006l )) still indicates 



that the best agreement is when there is an enrichment of heavy elements by between 0.1 
and 0.2 dex. 

We also conside red the high-ionizatio n color-color diagram because the [O III] 5007 A 



line was detected by 



Madsen et al. 



(120061 ) in several of their samples of Barnard's Loop. In 
each case there was significant disagreement with the models best fitting the low-ionization 
color-color diagram in the sense that there was excess [O III] emission. In the same way 
as M 43, this excess is almost certainly due to scattered light from the bright central H II 
regions in Orion where the I([0 III] 5007 A)/I(H/5) ratio is about three and a contamination 
of only a few percent is enough to account for the disagreement. 

As a guideline for future observations of Barnard's Loop, we present in appendix 
B the predicted relative intensities of the strongest emission lines over a wide range of 
wavelengths. We adopted T star = 31000 K, log $ = 7.91, and an abundance enrichment of 
0.15 dex as the best fit to the Barnard's Loop. 



3.6. Comparison of Barnard's Loop Properties with those of the 
Orion-Eridanus Bubble and the WIM 



Madsen et al. 



( 120061 ) measured a large number of samples within and near the 
Orion-Eridanus Bubble in addition to two samples (Wl and W2) that clearly lie within 
Barnard's Loop. They point out the continuity of the low-ionization color-color diagram 
results for all of these, in spite of a factor of 40 range of Ha surface brightness. Their results 
for the low surface brightness H II regions fall below the pattern established for the other 
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samples. In Figure 7 we show the data from Madse n et al.'s (2006 



Peimbert et al. 



igure 7 and have added 



(119751 ). The similarity of 



our new Barnard's Loop observations and those of '. 
the Barnard's Loop ratios (especially our own) and those in the Orion-Eridanus Bubble and 
the WIM samples is striking and probably indicates a common set of physical conditions 
and processes. 

To these observations we have added the predictions for our models used to explain 
the Barnard's Loop observations shown in Figure 5. The log $=7.9 calculations have 
been designated by the more general term, the ionization parameter U (log U=-3.07). 
U is the ratio of the number density of ionizing photons and hydrogen atoms and is 
given by U=0/(c n#), where c is the velocity of light and nn is the total hydrogen 
density. U is a broadly useful parameter because it directly indicates the conditions for 
photoionization, with various values of flux and gas density producing the same conditions 
for photoionization. U becomes a less useful measure in considering bright and dense 
photoionized gases because it does not give an indication of collisional processes. The same 
coding of log (Z/H) has been used for both Figure 5 and Figure 7 . In addition to the set 
of models displayed in Figure 5, we also show a series of calculations using our hottest stars 
at 40000 K (corresponding to a spectral type of about 07.5) at log U=-3.07 and stars of 
34500 K and 40000 K for log U=-3.67. 

Examination of Figure 7 indicates that variable abundance models for log U=-3.67 and 
T sta r=3100 K, log U=-3.67 and T star =40000 K, and log U=-3.07 and T star =31000 K fully 
enclose the space occupied by the Barnard's Loop, the Orion-Eridanus Bubble, and the 
WIM sample observations, with only the low surface brightness H II region observations 
falling below. The region occupied by the low surface brightness H II regions demands 
higher values of U than the Barnard's Loop although comparable values of Z/H and T sta r- 
Within the envelope of these bounding models one cannot tell from the low-ionization 
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color-color diagram of the WIM samples what causes an individual point to have its specific 
location, since the location is defined by U, T star , and Z/H. In the case of the Barnard's 
Loop samples, where the Z/H and T star values must be nearly constant, we'd expect 
variations in U to distribute the observed points along a line of constant color (a nearly 
vertical line) and that is the case. The wider distribution of the other Orion-Eridanus 
Bubble samples (the darker Madsen et al. (2006) points in Figure 7) would then indicate 
variations with position of abundance and/or T sta r of the dominant ionizing star, in addition 
to variations in U. Vari ations in T s t a r are certainly possible in such a large-scale sampling 



since there is evidence ( 10 'Dell 



20011 ) that the optically thick foreground Veil of the Orion 
Nebula is probably optically thin to the southwest and this would allow radiation from 
the hottest star in the region to illuminate Orion-Eridanus Bubble components in that 
direction. In the case of the WIM samples it is expected that there could be a significant 
range in photoionizing star temperatures, U, and possibly Z/H. 

If there is only a single value of the (Z/H) and it is that adopted for our M 43 
calculations then most of the Orion-Eridanus Bubble and WIM ratios can be explained 
by log U values between -3.07 and -3.67, with T star values of up to slightly more than 
35000 K. However, the lower left population of the Orion-Eridanus Bubble samples and the 
Barnard's Loop samples would require unrealistically low T star values, indicating that there 
must be regions of higher than average Z/H. 

We can constrain the likely T e of the Barnard Loop samples since they are all 
illuminated by the same radiation field. Table 2 gives the T e in the [N II] emitting zone 
for all our models. The two most closely matching the low-ionization color-color diagram 
are those with log U=-3.67 and log U=-3.07 with T star =31000 K, both with an abundance 
difference of 0.1 dex, and these have expected T e of 5970 K and 5940 K respectively. We 
will adopt a value of 5960±50 K for comparison with direct determinations. There is a great 
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uncertainty about the expected T c of the other parts of the diagram. For example, the 
T s tar=31000 K and log U=-3.07 model with average Z/H predicts line ratios about nearly 
the same as the T star =40000 K and log U=-3.67 model with an abundance enhancement of 
about 0.5 dex, would have the very different temperatures of 6530 K and 4680 K. 

The important conclusion of this section is that theoretically one can explain the 
low- ionization color-color diagram for the non-Barnard's Loop samples by a range of values 
of U, T star , and Z/H. A simple ratio of nebular line intensities for two different ions cannot 
produce an unambiguous estimate of T e . 



3.7. Direct Determinations of the Electron Temperatures in H II Regions, 

Barnard's Loop, and the WIM 

The reason for the designation WIM (warm ionized medium) is the fact that T e there 
is higher than the cold gas in the ISM. The actual value for the WIM's temperature is much 
more uncertain than sometimes stated in the literature because there are few observations 
(summarized in this section) that allow a direct determination and the indirect methods 
(based on only the I([N II] 6583 A)/I(Ha) ratio) commonly employed are uncertain as they 
assume a fixed nitrogen ionization ratio and abundance. In this section we summarize the 
results for T e derived by direct means. 

There are three direct methods of determining T c . The first is from the measurement 
of forbidden line intensity ratios within a single ion. The second is from the width of 
emission lines from ions of very different mass. The third is from the ratio of continuum to 
recombination line emission. 

Observations of the low-ionizationization WIM cannot use the most w i dely u sed T e 



indicator, the [O III] auroral/nebular line ratios. However, iReynolds et al.l (120011 ) were 
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able to measure the auroral 5755 A and nebular 6583 A lines of [N II] in a number of low 
surfac e brightness H II reg ions and several samples of WIM clouds along a single line of 



sight. 



Madsenetal. 



(120061 ) were able to measure these same line ratios in several additional 
WIM clouds and WIM clouds lying along the same line of sight as low surface brightness 
H II regions. They conclude that the WIM component has a line ratio corresponding to 
T about 2000 K higher than their sample of H II regions and that T is higher for WIM 
clouds of lower surface brightness. 



We show the iReynolds et al. 



f|200lh and 



Madsen et al. 



have added 



2005 



2006al ]bf) 



ine ratios for many additi onal b righter H II regi ons 



O'Dell k Harrisl (120101 ). The 



Reynolds et al. 



(120061) data in Figure 8 and 



Garcia- Roj as et al. 



(12004 . 



(1200 lh WIM sample had four 



velocity components and we see that when considered together (the "total" sample) the line 
ratios are similar or lower than most H II regions, indicating similar T e . Only the WIM 
component along the line of sight to Sivan 2 has a significantly higher inferred T e and it is 
comparable to the hottest classical H II region (NGC 3603). In contrast, the low surface 
brightness H II regions are systematically lower in line [N II] line ratio than the bright H II 



Reynolds et al. 



( 120011 ) sample have 



regions, hence have lower T e . The H II regions in the 
Ha surface brightnesses of 68-339 R. They are much lower surface brightness than classical 
bright H II regions. For reference, the O rion Nebula has an ex tinction corrected maximum 



surface brightness in Ha of 1.1 x 10 6 R (1 O'Dell & Harris 



2010h . 



he co nclusion that o ne ca n draw from the data presented in the 



( 120011 ) and 



Madsen et al 



Reynolds et al. 



( 120061 ) studies and comparison with the results from classical 



bright H II regions is that the only direct determinations of T e of WIM components 



indicate temperatures compar able to classica 



brightness H II regions in the IReynolds et al. 



aright H II 



(|200lh and 



regions. It is o nly t 



re low surface 



Mad 



are of unusually low T e , in spite of the emphatic statement in 



sen et al 



Haffner et al. 



(200 6J) sam ple that 



( 120091 ) that the 
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WIM temperatures are elevated to the bright classical H II regions. That statement is only 
true for the very low surface brightness components of the WIM, as we discuss in § 4.3. 
This can be due to them being photoionized by cooler stars. An additional factor may be 
that collisional de-excitation is not important at the lower densities in the WIM clouds so 
that cooling radiative transitions are relatively more important. There is an observational 
selection effect in the WHAM studies in that the low surface brightness H II regions are 
much larger than the typical H II regions since they are usually larger than the 1° diameter 
of the WHAM field of view. 

Line widths will characteristically have several components, the thermal width (which 
will scale as the square root of the ratio of temperature and ion mass) and random large 
scale mass motions along the line of sight. Since the atomic mass of S is 32 times that of 
H, one can hope to determine T e after making reasonable assumptions abo ut the common 



properties of the mass motion. In an early WHAM study, iReynoldsl ( 119851 ) compared the 
line widths of Ha and [S II] lines, but the spectral resolution and intrinsic widths of the 
lines did not allow an unambiguous determination of T e . Line widths have been measured 
in the diffuse H II region surrounding £ Oph, where the fainter components are comparable 
to the WIM in surface brightness. Observations of [N II] 6583 A, [S II] 6716 A, and Ha were 



interpreted by 



Baker et al. 



(120041 ) in a paper tha t only appears as an abstract. However, a 



Haffner et al 



(120091 ) . There one sees that 



summary of their results is shown as Figure 4 of 
most of the observed points lie between 6000 K and 9000 K and that there is a systematic 
and position dependent change in the non-thermal component. Again these temperatures 
are comparable to the classical bright H II regions. 

Comparison of optical emission lines and radio contin u um c an also provide a good 
source of T e . In the case of Barnard's Loop, iHeiles et al] ( I2OO0I ) use their own radio 
observations in two large samples at four wavelengths and WHAM optical surface 
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brightnesses to determine that T e is about 6100 K, in excellent agreement with the value 
of about 5960±50 K inferred from the low-ionization line ratios discussed in § 3.6. An 
unusual result using this approach occurs for the WIM 



Dobleretal. 



(120091 ) compared 



the continuum measured with the WMAP satellite with the Ha surface brightness. They 
determined that the observed ratio required a gas temperature of 3000 K. It is difficult 
to understand such a low temperature arising from photoionization processes, which 
has lead to the creation of a three component model composed of photoionized gas. 



gas that i s recombining and c oolin g, and cool neutra 



Following 



Wood & Revnoldsl (I1999h 



hyd 



rogen ( IDong &: Draine 



201lh . 



Dong fe Draind ( 1201 ll ) assumed in their model that 



there was a significant scattered light component, and found that a 15% scattered light 
contribution to Ha was necessary to produce a satisfactory model. It is hard to assess this 
multi-free-parameter model since it assumes that the warm gas component temperature 
is that indicated by the I([N II] 6583 A)/I(Ha) ratio method. There are also arguments 
based on direct obseryation s that high latitude clouds are strongly affected by scattered Ha 



radiation ( iWitt et al. 



in the detailed models of 



20101). If contamina t ion b y scattered Ha was stronger than derived 



Dong fc Draind (120111 ). then the temperature of the WIM gas 



contributing to the radio continuum would be higher than 3000 K. 



4. Discussion 



In this section we consider the results derived from the earlier sections. We consider 
the effects of variations of Z /H on the low-ionization color-color diagram, the general utility 
of using the low-ionization color-color diagram for determining T e , the systematic changes 
of T e determined from the low-ionization color-color diagram, and the issues involved with 
using large spatial scale line ratio variations from observations of other galaxies. 
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4.1. Effects of variations of Z/H on the low-ionization color-color diagram 

Our most important conclusion is that one cannot explain the low-ionization color-color 
diagram for the Barnard's Loop samples by the stars capable of causing this large-scale 
photoionization if the Z/H ratio is our reference value used in the M 43 calculations. We 
invoke a heavy element enhancement of about a factor of 1.4 to reconcile our models with 
the observations. A similar enhancement is also necessary to explain the high surface 
brightness population of the Orion-Eridanus Bubble samples. If such variations are 
necessary to explain the best studied samples, it is likely that the abundance variations also 
occur in the clouds producing the components of the WIM population. 

There certainly is evidence for local variations in Z/H within the ionized components 



of the interstellar medium. Two stu dies of multiple H I 



at various galactocentric distances (ISimpson et al. 



1995 



regions and compact 



Afflerbach et al. 



T II regions 



19971) found a 



general decrease in Z/H with increasing distance from the center of the Galaxy. Their 
analysis relies on ratios of infrared emission lines and will not be affected by uncertainties 
in the gas temperature or density. More important to our problem, they found in the local 
part of the Galaxy, variations of 0.3 dex both above and below the average and well beyond 
their estimated probable errors. These variations probably arise because stars with a range 
of masses and lifetimes produce different elements that enhance their locale. 

As noted in § 3.5, we have assumed in our models that the relative abundance of 
various elements vary together. For a small-scale object, where the products of an individual 
evolved star can be important, this assumption would be questionable. However, for the 
large objects that we consider here (Barnard's Loop, the Orion-Eridanus Bubble and the 
WIM clouds) the abundances must be determined by the products of multiple stars and the 
adoption of a constant ratio of individual elements is justified. 
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4.2. Application of the low-ionization color-color diagram for determining T c 



During the last decade there have been multiple papers on the WIM that argue 
that a low-ionization color-color diagram similar to Figure 7 can b e expl ained primarily 



by variations in T c as noted in t he recent revi e w by 



Haffner et al 



goes back to t 



he W IM study by 



( iMadsen et al. 



Haffner et al. 



fl2009j). The method 



(Il999l ) and has frequently been employed 



20061 ) . Within the assumption that the ionization ratio of nitrogen remains 
constant and a known N/H abundance applies, it is a simple matter to draw vertical lines 
representing constant values of T e in Figure 7. The nitrogen ionization ratio is defined 
as (H/H + )(N + /N), where H and N represent the total number density of hydrogen and 
nitrogen atoms and the superscripted values their ion number density This reflects the 
fact that the emissivity of a recombination line and forbidden line have reverse-sense 
dependencies on T e . However, if the nitrogen ionization ratio is lower than assumed the 
relative emissivity per unit volume decreases and a vertical line indicating a fixed T e would 
move to the left. If the N/H ratio is lower than assumed, the displacement would also be to 
the left. Of course both assumptions could be incorrect in different senses and one error can 
correct for the other, but without a good knowledge of both the nitrogen ionization ratio 
and the relative abundance of nitrogen and hydrogen, the method is suspect. 

In an attempt to identify the range of probable values of the nitrogen ionization ratio, 
we have extracted this information from our calculated models. For the models that most 
closely match the distribution of the Barnard's Loop, Orion-Eridanus Bubble, and WIM 
observations the ionization ratio varies little. The nitrogen ionization ratio for log U=-3.67 
and T sta r=31000 varies only from 1.007 to 1.014 over the range of Z/H from -0.5 to 0.5 
dex. The ionization ratio for log U=-3.07 and T sta r=40000 varies only from 1.012 to 1.035 
over the range of Z/H from -0.5 to 0.5 dex. This confirms that variations in the ionization 
ratio do not play an important role, as previously assumed and calculated. This is in 
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excellent agreement wi t h the predictions of Sembach et al. (2000). Unfortunately, in the 



study of 



Madsen et al 



( 120061 ). which drew on the Sembach et al. (2000) models a value 
of the nitrogen ionization ratio of 0.8 was adopted, whereas this is actually the value for 
N + /N. This error w as not corrected when the results were repeated in a review article 



( jHaffner et al 



20091). This makes the electron tempera tures they pre s ent to be too large 



Haffner et al. 



(119991 ) a nitrogen 



for their assumed abundance. In the original study of 
ionization ratio of 1.0 was adopted, which means that those temperatures should be correct, 
if the nitrogen abundance they adopted of N/H=7.5 x 10 -5 is both correct and uniform. 

The more important limitation of the I([N II])/I(Ha) ratio method is the other scaling 
factor, the N/H ratio. We argue in § 3.6 that in the case of Barnard's Loop that there is a 
Z/H abundance enhancement of about 0.15 dex that (alone) would shift a line of fixed T e 
to the right by a fac t or of 1.4 in Figure 7. As noted above, variations of this magnitude 



(jSimpson et al 



1995 



Amerbach et al 



19971 ) are known to exist. Without a knowledge of 



the abundance, which is usually derived for gaseous nebulae after one knows T e by direct 
means, it is impossible to determine an accurate value of T e for specific objects from 
low-ionization color-color diagrams. 



4.3. Systematic variations of T c determined from the low-ionization 

color-color diagram 

Because the progression of calculated points for abundance excesses or deficits makes 
a loop that peaks near the predictions for the nominal abundance, abundance variations 
are unlikely to explain the high values of the low-ionization ratios of the lowest surface 
brightness WIM components. These components must have higher values of T e . These 
high value components can be explained by photoionization processes if the illuminating 
stars are of higher T star or are stars of lower temperature whose LyC radiation field mimics 
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that of a hotter star (IWood fe Mathial2004l ). whereas iReynolds et al.l ( 119991 ) argued that 
photoionization processes are insufficient. Modifying the photoionizing radiation is possible 
through selective removal of photons of energies slightly greater than the ionization energy 
of hydrogen of 13.6 eV, a process commonly known as radiation hardening. If one introduces 
non-photoionization processes that come into play at these very low densities, then one 
doesn't need to assume a harder radiation field. 

In summary we can say that the low-ionization color-color plots of Barnard's Loop, 
the Orion-Eridanus Bubble and the higher surface brightness components of the WIM can 
be explained by combinations of T e and abundance variations, using available stars. In the 
case of the lowest surface brightness components of the WIM it is necessary to either modify 
the characteristic radiation field or to introduce non-photoionization heating processes. 



4.4. Large scale line ratio variations in other galaxies 



Study of the diffuse ionized gas in other galaxies may provide some help in 
understanding our own WIM and the cause of the systematically higher values of T e in 
the lowest surface brightness components. In the study of other galaxies one has the 



advantage of easily loo 



(jTullmann &: Dettmar 



ring fo r yariations with position and these are commonly found 



2000a 



b; 



Otte et al 



2001 



20021 ) . but this is at the expense of 



losing the diagnostically valuable tool of being able to divide the contributors into surface 
brightness groups and it often appears necessary to invoke non-photoionization processes 
to explain the observations. The work that we report on here builds from the physics 
operating in a succession of photoionized objects of decreasing density and increasing scale 
(M 43, Barnard's Loop, the Orion-Eridanus Bubble, and local components of the WIM. 
These final conclusions can then be a jumping-off point in discussion of the diffuse ionized 
gas in other galaxies. 
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5. Conclusions 

We have been able to reach several important conclusions from this study that began 
with new observations of the Barnard's Loop. These have provided data similar to previous 
studies but in a new region of the object and at higher spectrophotometric accuracy. 
These observations were supplemented by intensive photoionization modelling. The major 
conclusions are: 

1. Barnard's Loop is photoionized by the most luminous stars in the Orion constellation 
except for 9 l Ori C and # 2 Ori A, whose radiation is largely absorbed locally. 

2. Tests of our photoionization models on the recently well observed and intrinsically 
simply H II region M 43 give a good fit to the low-ionizationization primary line ratios for 
stellar models close to the T sta r of the complex exciting star NU Ori. 

3. Barnard's Loop is similar in its properties to other regions in the Orion-Eridanus 
Bubble and lies at the high surface brightness end of a population of components of the 
WIM. 

4. Our best models that explain the low-ionizationization lines fail to predict the 
observed strength of the [O III] lines in both M 43 and the Barnard's Loop. This is due to 
a small contamination of the nebular emission by scattered light arising from M 42. 

5. The location of the Barnard's Loop observations in a low-ionization color-color 
diagram cannot be explained from photoionization by the most likely dominant stars unless 
one assumes a local Z/H enhancement of about 0.15 dex. This argues that Barnard's Loop 
is enriched and possibly shaped by high velocity mass loss from evolved stars near its center. 

6. The electron temperature derived from our models of Barnard's Loop and the 
low-ionization line ratios i s about 5960±50 K, in excellent agreement with the optical/radio 



method result of 6100 K (IHeiles et al 



2000) 
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7. The population of low-ionization color-color line ratios of Orion-Eridanus Bubble 
and WIM components is enclosed by a small range of values of U, T star , and Z/H. 

8. Comparison of T e derived by direct methods in classical high surface brightness 
H II regions and a few samples of the WIM indicates that the WIM components are 
of comparable T c to H II regions . Only very low surface brightness H II regions have 
systematically different and lower T c . 

9. The lowest surface brightness components of the WIM in the low-ionization 
color-color diagram are likely to have systematically higher T c than the higher surface 
brightness components. 

10. We establish that the usual method of determining T e in WIM components through 
I([N II])/I(Ha) ratios is subject to an important uncertainty arising from known abundance 
variations. 

11. Small changes of abundance from the reference value produce a confined short loop 
in the low-ionization color-color diagram for a fixed T star . This means that if there is only 
a variation in abundance within ±0.3 dex, the lowest surface brightness components of the 
WIM are systematically illuminated by harder radiation fields and that photoionization 
processes can explain the observed line ratios. 

Partial financial support for GJF's work on this project was provided by National 
Science Foundation grants AST 0908877 and AST 0607028 and National Aeronautics and 
Space Administration grant 07-ATFP07-0124. CRO's work was partially supported by 
STScI grant GO 10967. 

Facilities: CTIO(1.5 m). 
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A. Recent Changes to Cloudy 



This paper makes extensive use of two newly-introduced facilities. Cloudy can now 



interpolate upon large grids of stellar atmosphere spec tral energy distri 



of the TLUSTY OSTAR2002 and BSTAR2006 SEDs flLanz fc Hubenv 



Dutions. In the case 



2003 . 



20071 ) we can 



interpolate on effective temperature, surface gravity, and metallicity, as we cannot vary 
abundances individually. We used TLUSTY to be con sistent with the physical calibration 



of the spectral classes presented by 



Heap et al. 



(120061 ). The interpolation methods have 



been generalized and many other grids of stellar SEDs are available. We have developed 
a domain decomposition method to compute these grids on MPI-aware parallel machines. 
Each grid point is an independent model calculation and so can be done on separate 
computer nodes. This results in a speedup that is of the order of the number of available 
processors. 

The capability to compute grids of photoionization models where certain key parameters 



are incremented in equidis tant linear or 



and has been discussed in 



Porter et al 



ogarithmic steps was introduced several years ago 



(120061 ). We recently enhanced this capability by 



parallelizing the algorithm using the Message Passing Interface (MPI) specification. This 
allows much larger grids to be computed in parallel on distributed clusters of computers. 
The calculations are set up in such a way that each core calculates a separate model 
(domain decomposition) and all the results are gathered when the grid has finished. The 
communication overhead is negligible since the MPI threads only need to communicate 
when the grid calculation starts and finishes. Hence this algorithm is highly efficient and 
scales well to high numbers of cores for sufficiently large grids. 



B. Predicted Line Ratios 
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Table 1: Flux Rat 


ios for the Barnard's 


Loop Samples. 






Wavelength (A) 


Source 


This Study 


Peimbert et al. 1975 


WHAM1 


WHAM2 


3727 


[0 II] 




0.60 






4340 


H 7 


0.51 


0.51 






4861 


H/3 


1.00 


1.00 






5007 


[0 III] 


<0.007 


<0.22 


0.020 


0.004 


5876 


Hel 


<0.007 




0.015 


0.004 


6548 


[N II] 


0.20 








6563 


Ha 


2.90 


2.66 


1.00 


1.00 


6583 


[N II] 


0.66 


0.68 


0.23 


0.23 


6716 


[S II] 






0.13 


0.16 


6716+6731 


[S II] 


1.00 


0.79 
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Table 2: Predicted Electron Temperatures in the [N II] Emitting Zone. 




Log U=-3.67 


log U=-3.67 


log U=-3.07 


log U=-3.07 


A log (Z/H) 


T star =31000 K 


T star =40000 K 


T star =31000 K 


T star =40000 K 


-0.5 


8950 


10280 


9270 


10940 


-0.4 


8490 


9890 


8730 


10420 


-0.3 


8020 


9440 


8200 


9860 


-0.2 


7530 


8950 


7660 


9270 


-0.1 


7010 


8410 


7100 


8650 


0.0 


6500 


7850 


6530 


8020 


0.1 


5970 


7260 


5940 


7360 


0.2 


5420 


6670 


5330 


6680 


0.3 


4850 


6040 


4700 


5970 


0.4 


4210 


5380 


4000 


5240 


0.5 


3480 


4680 


3280 


4440 
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Table 3: 


Predicted Line Ratios. 


Ion 


Wavelength 


T /TYW/^ 
L/L{tlp) 


Lya 


IzlO A 


lZ.4Do 


C II] 


2620 A 


U.ol4o 


Mg II 


970.8 A 
2 I »o A 


n 07 AO 

u.z i 4y 


[0 II] 


o/zV A 


o 1 no 


[S II] 


/i n7n A 
4U/U A 


O OQSO 

U.Uooy 


[S II] 


AO 78 A 

4U 1 O A 


n m 
u.uiou 


H/3 


4o61 A 


l.UUUU 


[N II] 


l DO A 


n m no 
u.uiuy 


He I 


c;a7R A 

Do l A 


n m qi 
u.uiy i 


[0 1] 


conn A 
OoUU A 


U.UlOo 


[N II] 


0D4o A 


U.Z14o 


Ha 


OOOO A 


z.oyuo 


[N II] 


fi^8/1 A 
0004 A 


U.OOO ( 


[S II] 


fi71 £ A 
( 10 A 


U.OoOl 


[S II] 


R7Q1 A 


ri 9707 
0.2 ( UV 


[0 II] 


7QOQ A 

1 626 A 


U.uz4y 


[0 II] 


7QQO A 


n n,9n f ? 
u.uzuo 


[S III] 


QflfiQ A 
yuoy a 


n n7^ 

U.U ( JO 


[S III] 


0^*39 A 
yooz a 


n 1 8RQ 
u. looy 


He I 


1 08^ //m 
l.UOO /Alii 


u.uzyo 


H I 


1817 //m 


n m 9^ 


[Ar II] 


fi Q80 //m 

u.you /im 


u.uouu 


[Ne II] 

L J 


12.81 pm 


0.1425 


[S III] 


18.67 /im 


0.0824 


[S III] 


33.47 /im 


0.1916 


[Si II] 


34.81 /im 


0.0202 


[N II] 


121.7 fan 


0.0306 


rc in 


157.6 um 


0.3029 
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Fig. 1.— This 24?06 x 16?06 groundbased telescope image with north at the top has 
superi mposed the posit i ons o f the slit used for our spectroscopic study(CTIO-2008) and 



that of 



Peimbert et al. 



the study of 



(119751) (P RT-P R7a), in addition to the two WHAM apertures in 



Madsen et al. 



fl2006D(WHAM 1 and WHAM2). The circle around the CTIO 
slit positions are there only to aid finding their locations since this long (for a spectro- 
graph) slit is small when looking at an object as large as Barnard's Loop. The stars 
most luminous in Lyman Continuum (LyC) radiation are identified . The line indicating 



10 pc length assumes a distance of 440 pc ( 10 'Dell &: Henney 



20081 ) . The filter used for 



the image isolated the Ha+[N II] lines and is used with the permission of Peter Erdmann 



(http:/ /messier.obspm.fr/xtra/ngc/b-loop.html). 
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Fig. 2. — The final spectrum of our new sample of Barnard's Loop is plotted in arbitrary 
units. Only an upper limit of about 0.007 for the ratio of flux relative to the H/3 line can 
be established for the [O III] 5007 A. The negative signals are the result of over-subtracting 
the sky [O I] lines, as discussed in the text. 
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Fig. 3. — These two color-color plots present the results of our calculations for log 0(H) = 
10.4, 10.7, and 11.0 (photons cm" 2 s _1 ), appropriate for M 43, with a range of T star of 27542 
K through 40040 K in steps of 0.0325 dex, as described in th e text. The open circle s with 



O'Dell k Harrisl (120101 ). The 



error bars give the results for the six M 43 sample spectra in 
open squares surround the point calculated for T star =29682 K, the closest sample to that 
expected from the B0.5V spectral type of NU Ori. Using the known approximate effective 
temperature of the star, we see that the most likely solution is log</>(H) = 11.0. 
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Fig. 4. — The M 43 predicted T e for the stellar atmosphere models and the best fitting flux 
(log0(H) = 11.0) from Figure 3 are presented. Filled squares are the predictions for [S II] 
and the filled circles are for [N II] . The average T e values and their dispersions for M 43 as 



derived from the auroral to nebular line ratios of [N II] and [S II] by 
are also shown. 



O'Dell & Harrisl (120 101 ) 
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Fig. 5. — This low-ionization color-color plot is similar to Figure 3 (left). The series of 
open circles is a set of calculations using the lowest T star in our calculations (27500) for 
various values of log $ ( 6.7, 7.3, 7.6, 7.9, 8.2, 8.5, and 9.1 top to bottom, with the 7.9 
value, favored by the model of the Barnard's Loop, shown in grey). The progression of filled 
circles are predicted line ratios for a series of heavy element to hydrogen abundance ratios 
(Z/H) in steps of log (Z/H) =0.1 relative to the values employed in the M 43 calculations 
(enclosed by an open square). The extreme values of the log (Z/H) differences are labeled. 
For this sequence of calculations we have used the well defined parameters (T star = 31000 
K, log $=7.9) for the Barnard's Loop. We have added the observed results for this study 
(BL), Peimbert, et al.'s 1975 study (P7 5), and the two WHAM points (Wl and W2) lying 



on Barnard's Loop (IMadsen et al. 



20061 ). In each case the unc ertainty of the observ ed point 



Madsen et al. 



(J2006J) data by 



is within the size of the enclosing ellipse. A correction of the 
a factor of 1.67 has been made since the WHAM [S II] data are for the 6716 A line only. 
The progression from low Z/H to high Z/H is a progression towards lower T e . The results 
shown as open circles makes it clear that no values of log $ with log (Z/H) = matches the 
constraints imposed by the observed line ratios. For ease of comparison, the range of values 
shown are the same as in Figure 7. 
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Fig. 6. — This high-ionization color-color plot is similar to Figure 3 (right) in the selection 
of the observed line ratios. The progression of predictions for varied values of log(Z/H) in 
steps of 0.1 dex with log 0(H)=7.9 (log U= -3.07) is shown as a series of filled circles with 
connecting lines. The largest log(Z/H) difference is labeled and the square encloses the value 
for the average Z/H model. The Barnard's Loop observations are shown with the same 
symbols as in Figure 5. In each case the uncertaint y of the observed po int is within the 



Madsen et al. 



(2006) observations of 



Madsen et al. 



fl2006h 



size of the enclosing ellipse. The open circles are the 
regions in the Orion-Eridanus Bubble and the filled squares are the 
observations of low surface brightness H II regions. In these cases the errors are less certain 
and probably not as great as the dispersion of the value of the samples. Upper limits are 
indicated with a descending arrow. The region occupied by the low surface brightness H II 
region observations generally correspond to higher values of the ionization parameter than 
used in the Barnard's Loop calculations. 
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Fig. 7. — This figure is like Figure 5 but with additional obse rvations of WIM components 



and additional theoretical models. The 



Madsenetal 



( 120061 ) observations of low surface 



brightness H II regions are shown as filled squares. The data from Madsen et al.'s (2006) 
Figure 7 depicting ratios in the vicinity of the Orion- Eridanus Bubble are shown with error 
bars. The darker Madsen points to the left are higher surface brightness samples closer to the 
Galactic plane and falling within the Orion- Eridanus Bubble, while the lighter Madsen points 
are lower surface brightness samples further from the Galactic plane and are considered 
samples of the WIM. We also show the predictions for our photoionization models with 
variable abundances for stellar temperatures of 31000 K and 40000 K for two values of 
the ionization parameter (the upper pair of calculations for 31000 K and 40000 K are for 
log U=-3.67 and the lower pair of calculations for 31000 K and 40000 K are for log U=-3.07). 
The connecting lines for 31000 K and log U=-3.07 vary in color for clarity. The pattern of 
variation of differences in log (Z/H) within a sequence is the same as in Figure 5. 
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Fig. 8. — This figure presents the only published direct measurements of T e dependent [N II] 
line ratio in the WIM and compares it with the same line ratio in well-known H II regions. 
Larger y values correspond t o higher T P . Open diam onds are the individual velocity co mpo- 



nents of the WIM samples of 



Reynolds et al. 



(l200lh (C1-C4) and 



Madsen et al. 



(120061 ) (Arc 



and WIM(Sivan 2)), the filled diamond is the result of considering all velocity components 
of the C1-C4 samples as a single line, and the open squares are WHAM measurements of 
l ow surface brightness H II regions. Filled circles are for the Ori on objects M 42 and M 43 



O'Dell k Harris 



2005 



2006a 



20101 ). and open circles are from the studies of iGarcia-Rojas et al. 



(12004 



bj). There is no indication that the WIM components are systematically higher 
T e than the higher density well-known H II regions. The dashed lines on the right indicate 
the I([N II] 5755 A)/I([N II] 6583 A) ratios expected at various T e . 



